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Abstract 

Since the smallest leptonic mixing angle 6 13 has been measured to be relatively large, it is quite 
promising to constrain or determine the leptonic Dirac CP-violating phase 5 in future neutrino 
oscillation experiments. Given some typical values of 5 = tt/2, tt, and 3n/2 at the low-energy scale, 
as well as current experimental results of the other neutrino parameters, we perform a systematic 
study of the radiative corrections to 5 by using the one-loop renormalization group equations in 
the minimal supersymmetric standard model and the universal extra-dimensional model. It turns 
out that 5 is rather stable against radiative corrections in both models, except for the minimal 
supersymmetric standard model with a very large value of tan /3. Both cases of Majorana and 
Dirac neutrinos are discussed. In addition, we use the preliminary indication of 5 = (1.08^0 31) 7r 
or 5 = (1. 67jlo 77) n horn the latest global-fit analyses of data from neutrino oscillation experiments 
to illustrate how it will be modified by radiative corrections. 
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I. INTRODUCTION 



In the last two decades, our knowledge on neutrinos has been greatly improved by a 
number of elegant neutrino oscillation experiments Now, we are convinced that neutrinos 
are massive, and they can transform from one flavor to another when propagating in vacuum 
or in matter. The lepton flavor mixing phenomenon can be described by a 3 x 3 unitary 
matrix V, namely the leptonic mixing matrix, which is conventionally parametrized through 
three mixing angles 9 l2 , 9 l3 and 9 23 , as well as three CP- violating phases 5, p and a, viz., 
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where = sin 9^ and = cos 9^ for ij = 12,13,23. Note that P = diag(e ip , e 1<T , 1) is 
a diagonal matrix with p and o being two Majorana-type CP-violating phases if neutrinos 
are Majorana particles, while P = 1 if neutrinos are Dirac particles. Current experimental 
data indicate that the three leptonic mixing angles are 9 12 ps 34°, 9 13 w 9° and 9 23 ~ 40°. 



Two independent neutrino mass-squared differences are found to be Am 



7717, 



ml PS 



7.5 x 10" 5 eV 2 and |Am§! 



\m\ — m\\ ps 2.5 x 10~ 3 eV 2 . The latest global-fit results of 
neutrino parameters are shown in Table I. However, we are still unclear whether the neutrino 
mass ordering is normal (i.e., Am^ > 0) or inverted (i.e., Am^ < 0), and the leptonic Dirac 
CP-violating phase 5 remains experimentally unconstrained. 

The recent results from Daya Bay jsj and RENO jsj reactor neutrino experiments have 
established that 9 13 ps 9°, which is rather large. Hence, it is quite promising to determine the 
leptonic Dirac CP-violating phase 5 by comparing the oscillation probabilities of neutrinos 
and antineutrinos in future long-baseline neutrino oscillation experiments 4J. In addition, 
the km 3 -scale neutrino telescopes (e.g., IceCube and KM3NeT) could provide us with useful 
and complementary information about leptonic CP violation by precisely measuring the 
flavor composition of ultrahigh-energy astrophysical neutrinos |5j. If the Deep Core of the 
IceCube detector is made denser to lower the energy threshold down to a few GeV, such as 
the proposal PINGU [6|, a large amount of atmospheric neutrino events can be collected and 
used to determine the neutrino mass hierarchy and perhaps the leptonic CP-violating phase 
[3] • On the other hand, a lot of neutrino mass models based on discrete flavor symmetries or 



TABLE I: The best-fit values and la ranges of the neutrino parameters from the latest global-fit 
analyses of neutrino oscillation experiments, where the normal neutrino mass hierarchy is assumed. 
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Ref. [13] 
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sin 0\2 


0.291 - 0.325 


0.287-0.313 


0.303 - 0.336 




0.0241 
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0.0246 


sin 2 


0.0216- 0.0266 


0.0207- 0.0253 


0.0218 - 0.0275 




0.386 


0.410 


0.427 


sin 2 #23 


0.365-0.410 


0.385 - 0.447 


0.400 - 0.461 


Am^/ICT 5 eV 2 


7.54 
7.32 - 7.80 


7.50 
7.32-7.69 


7.62 
7.43-7.81 


Am^/ICT 3 eV 2 


2.51 
2.41 - 2.57 


2.47 
2.40-2.54 


2.55 
2.46-2.61 


8/tt 


1.08 
0.77-1.36 


1.67 
0.90-2.03 


0.8 
0-2.0 



phenomenological assumptions have recently been proposed to describe the observed leptonic 
mixing pattern, in particular a relatively large 8 13 . Interestingly, the leptonic CP- violating 
phase o has been predicted in some models to be rather large (e.g., 5 > tt/3) or even maximal 



i.e. 



S 



9[. In other models, leptonic CP violation is shown to be absent, namely 



le latest global-fit analyses of neutrino 



ll| and 5 



121 ]. although 



5 = or Ti [10( 1 . It is worthwhile to mention that t 
oscillation experiments yield 5 = (1.08l2;|f) 
the la errors are still quite large. 1 Therefore, we have already obtained some preliminary 
information on the leptonic CP-violating phase 5 from the global-fit analyses. 

In this work, we are concerned with how the theoretical predictions or the observed value 
of S will be modified by the radiative corrections when running from a low-energy scale to 
a superhigh-energy scale. This question does make sense if we believe that there exists at 
some superhigh-energy scale a unified theory for flavor mixing and CP violation in both 
quark and lepton sectors. Once the leptonic CP-violating phase 5 is measured in future 
neutrino oscillation experiments, the renormalization group (RG) evolution of 5 will tell 
us how large or small it will be at a given superhigh-energy scale. As a matter of fact, 



1 The best-fit value is found to be 6 = 0.8 ir for the normal mass hierarchy and 5 = —0.03 ir for the inverted 
mass hierarchy by another global- fit group [l3| . However, there is no constraint on 5 within the la range. 
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;he running of leptonic mixing parameters has been extensively discussed in the literature 
3 [ill, an d more recently in Ref. 16], where the authors concentrate on the newly measured 
8 l3 . Different from the previous works, we focus on 5 and perform a systematic study of 
its running behavior in the minimal supersymmetric standard model (MSSM) and in the 
universal extra-dimensional model (UEDM). The motivation for such a study is two- fold: 
(1) The leptonic CP-violating phase 5 is the last fundamental parameter (except for the 
neutrino mass hierarchy) to be measured in the future neutrino oscillation experiments, and 
now both the theoretical models and the global-fit analysis can provide us with preferred 
values of 5 at the low-energy scale. (2) The models with supersymmetry or extra spatial 
dimensions are the most natural extensions of the SM, which can solve the gauge hierarchy 
problem and offer good candidates for the dark matter. 

In lack of a complete theory for neutrino mass generation, we implement the dimension- 
five Weinberg operator to account for tiny Majorana neutrino masses 17| . The RG running 
of S in the case of Dirac neutrinos will be considered as well for comparison and completeness. 

The remaining part of the present paper is organized as follows. In Sec. II, we set up the 
basic framework for the RG running of leptonic mixing parameters in the case of Majorana 
neutrinos. The renormalization group equation (RGE) of S is derived analytically, and solved 
numerically. Section III is devoted to the RG running of 5 in the case of Dirac neutrinos 
in the MSSM. We summarize our conclusions in Sec. IV. The complete set of RGE's in the 
SM, MSSM, and UEDM for Majorana neutrinos are collected in Appendix A, while those 
in the SM and MSSM for Dirac neutrinos in Appendix B. 



II. RUNNING OF CP- VIOLATING PHASE: MAJORANA NEUTRINOS 

First of all, we derive the RGE for the leptonic CP-violating phase 5, assuming that 
neutrinos are Majorana particles. Without loss of generality, we introduce the dimension- 
five Weinberg operator responsible for neutrino masses 

-C u = ^(lH)-K-(H T £ c )+h.c, (2) 

where i and H stand for the lepton and Higgs doublet fields, respectively, and k is a symmet- 
ric and complex matrix of the inverse mass dimension. After electroweak symmetry breaking, 
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TABLE II: Explicit expressions of Re (U^U)^ 
parametrization of leptonic mixing matrix. 



and Im 



(WU)^ for i < j in the standard 



Re 



Im 



11 

22 
33 

12 9 12 + s 13 s s 9 2 3 

13 — Si2 c l3#23 + c 12 c 5^13 — C 12 S 13 C 13 S (5" 

23 +Ci 2 C 13 ^23 + s 12 c ,5^13 ~ S 12 S 13 C 13 S <5" 



+2s 12 c 12 s 13 s (5 6' 2 3 + c\ 2 s\ 3 8 



'12 l '12' : '13' : '<5 1/ 23 



-12 



2y ,s 13 s <5^23 "I" s 19. c 19. s 13^ 



12 u 12°13'- 



c 12 s <5^13 C 12 S 13 C 13 C <5^ 

— s 12 s (5 6 l 1 3 — s 12 s 13 c 13 c g o 



the mass matrix of three light Majorana neutrinos is given by M v = kv 2 with v ~ 174 GeV 
being the vacuum expectation value (vev) of the SM Higgs field, or by M v = k(v sin j3) 2 
with tan (3 being the ratio of the vev's of two Higgs doublets in the MSSM. Note that we are 
working within an effective theory, and consider the running of neutrino mixing parameters 
below the cutoff scale A where new physics takes effects. 
At one-loop level, the evolution of k is governed by 



IGvr 2 ^ = a K + C K 
dt 



(3) 



where t = ln(/i/A EW ) with \x being an arbitrary renormalization scale between the elec- 
troweak scale A EW ^ 100 GeV and a cutoff scale where new physics comes into play, and Y l 
is the Yukawa coupling matrix of the charged leptons. The coefficients a K and C K are flavor 
universal, and have been explicitly given in Appendix A for the SM, the MSSM, and the 
UEDM. It is worth stressing that Eq. (3) takes on the same form in all the models under 
consideration. However, the coefficients in the RGE's may differ. We will distinguish them 
by adding the corresponding superscripts to these coefficients, as shown in Appendix A. 

A. Analytical Results 

Since the RGE's of neutrino mass matrix M v = kv 2 in the SM and UEDM, or M v = 
k(v sin ft) 2 in the MSSM, are given by the same formula in Eq. (3), the evolution of neutrino 
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TABLE III: The coefficients IZfj and Xf- for a = e, /i, r and ij = 12, 13, 23 in the standard 
parametrization of leptonic mixing matrix. 
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mass eigenvalues and leptonic mixing parameters can be figured out in the same way. In 
flavor basis, where the Yukawa coupling matrix of the charged leptons is diagonal, namely 
Y { = D l = dia.g(y e ,y fl ,y T ), k can be diagonalized by the leptonic mixing matrix V, namely 
V^kV* = k = diafff k 1; k 2 , k, 3 ). Generally speaking, an arbitrary 3x3 unitary matrix V can 
be factorized as V = QUP, where Q = diag(e 1 ^ >e ', e 1 ^, e 1 ^) and P = diag(e ip , e 1<J , 1) are pure 
phase matrices, while the unitary matrix U consists of three mixing angles 6 12 , 9 13 , 6 23 and 
the Dirac CP- violating phase 5 [cf. Eq. (1)]. Although the phases <j) a (for a = e, /i,r) are 
unphysical and can be removed by rephasing the charged-lepton fields, we will keep them in 
the derivation of the RGE's for neutrino masses and leptonic mixing parameters. 

Since y\ < ^ < j/J, we take into account the dominant contribution from the tau-lepton 



Yukawa coupling to the RGE of k. Following Ref. 



, one obtains 



16^ = ^K + 2C^ Ti | 2 ) , (4) 

where a K and C K should bear the corresponding superscripts when Eq. (4) is applied to a 
specific model. Given m i = k^ 2 (for % = 1,2,3), we observe that Eq. (4) determines the 
evolution of absolute neutrino masses. Moreover, it is straightforward to find that U ai , p, a, 
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and 4> a (for a — e, p, r and i = 1,2, 3) have to fulfill the following equations: 
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where Qj = 4/t i K J /(/t 2 — k 2 ) and = 2(/t 2 + k 2 )/(k 2 — k 2 ) have been defined, and the 
overdot refers to the derivative with respect to the running parameter t. In addition, IZfj = 
Re {UaiU a j) and Z-* = Im (U* i U a j). Given the standard parametrization of U in Eq. (1), 
the matrix elements of UW are shown in Table II, while the coefficients 1Z? 4 and Xf- are 
given in Table III. Note that Eqs. (5) and (6) form an array of differential equations linear 
in {0 12 , 9 13 , #23, 5, p, a, e , , r }, which can be explicitly solved. As a result, the RGE of 5 
can be approximately written as 
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Since the last two terms in the third line of Eq. (7) are proportional to s 13 and s"{ 3 , we have 
neglected the terms further suppressed by (^(Am^/IAm^l). If neutrino masses are nearly 
degenerate m 2 ^> lAmU ^> Am^, which will always be assumed in the following, we have 
dj ~ dj ~ 4m ! 2 /(m^ — m 2 ) and C21 ^ IC32MC31I ^ 1; an d thus, the RG evolution of 5 could 
be significant. To next-to-leading order, Eq. (7) approximates to 
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where we have taken m l as the absolute neutrino mass and ignored the difference between 
Amgj and Am| 2 - Some comments are in order: 



In general, the evolution of 5 is dominated by the leading-order term — C2i s 23 s 
the right-hand side of Eq. (7). At higher order, if the terms suppressed by IC31I/C21 



2(p-a) 



on 
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1/30 are taken into account, then those by s\ 3 » 1/40 should also be 



kept for consistency, since they are of the same order of magnitude, as we have done 



in Eq. (8). The relative error in Eq. (7) is at the level of S13IC31I/C 
the best-fit values of 6 l3 and neutrino mass-squared differences. 



0.5 %, 



given 



It is evident from Eq. (7) that the evolution of 5 is entangled with that of three mixing 
angles and two Majorana CP-violating phases. In particular, it depends crucially on 
the Majorana phases p and a. It has been found that the Dirac CP-violating phase S 
can be radiatively generated from p and a, even if the initial value of 5 is vanishing 



19[. On the other hand, the RG evolution of S becomes negligible when p ~ cr, while 



the mixing angle 8 12 is quite sensitive to the RG effect in this case. 

• The RGE's of 5 in the SM, the MSSM, and the UEDM are given by the same formula 
in Eq. (7), but with different values of the coefficient C K . We have C^ M = —3/2 in 
the SM, while C^ SSM = 1 in the MSSM and C^™ = -3(1 + s)/2 in the UEDM, 
respectively. Therefore, given the same Majorana CP-violating phases and leptonic 
mixing angles, the evolution of 5 in the MSSM will be in the direction opposite to that 
in the SM and the UEDM. 

Finally, we observe from Eq. (5) that the identity <fi e + <fi + <fi T + p + & = holds in 
the standard parametrization of U. The proof is as follows. Given a general non-singular 
matrix X, whose elements are functions of the running parameter t, one can prove that 
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d[det(X)]/dt = det(X) • tr[X _1 (dX/di)]. If we take X to be a unitary matrix U with 
det(U) = 1 and U^ 1 = U\ then tr(WU) = can be obtained. This observation together 
with Eq. (5) leads to the identity e + M + <p T + p + & = 0. However, this identity depends 
on the specific parametrization of U. For instance, if det(£7) = with being the Dirac 



CP- violating phase, then we have <p = <p e + <p +<p T + p + &, as shown in Ref. [18]. 



B. Numerical Results 

We proceed in this subsection with the numerical solution to the RGE of the leptonic 
Dirac CP-violating phase 5. Since the evolution of S in the SM is negligible even in the case 
of a nearly-degenerate neutrino mass spectrum, we consider only the MSSM and the UEDM. 
Note that no approximations to the RGE of 5 will be made in our numerical calculations. 
Our numerical results are shown in Fig. 1, and the main points are summarized as follows. 

In the MSSM, we have taken two typical values of tan (3 = 10 and tan /3 = 30 for 
illustration. In both cases, the absolute neutrino mass m 1 = 0.1 eV is assumed, which 
is consistent with the cosmological bound m 1 + m 2 + m 3 < 1.3 eV (95 % C.L.) from the 



WMAP Collaboration 20|. For the initial values of 5 at the electroweak scale, we have chosen 
5 = it/2, it, and 3tt/2 as typical examples. Since the tau-lepton Yukawa coupling is given by 
y 2 T = rn? T {\ + tan 2 fi)/v 2 in the MSSM, the evolution of 5 should be significantly enhanced for 
a large value of tan (3, as shown in the upper plots of Fig. 1. For tan j3 = 30, the RG running 
of S is quite significant. In particular, even if § = tt is found at the low-energy scale, namely, 
there is no CP- violating effect in neutrino oscillation experiments, the maximal CP- violating 
phase 5 = tt/2 or Zir/2 can be achieved at the cutoff scale A = 10 14 GeV. In other words, 
one can change from the scenario with a zero CP-violating phase to that with a maximal 
CP-violating phase, or vice versa. For tan/3 = 10, the radiative correction to 5 is at most 
10 % even at A = 10 14 GeV. 

In the UEDM, we have input two different values of the absolute neutrino mass m l = 
0.1 eV and m 1 = 0.5 eV. As shown in the lower plots of Fig. 1, 5 is rather stable against 
radiative corrections for m 1 = 0.1 eV. Even for m l = 0.5 eV, which is marginally in tension 
with the cosmological bound, the relative change of S at the cutoff scale A = 3 x 10 4 GeV 
is not larger than 10 %. The cutoff scale A = 3 x 10 4 GeV in the UEDM has been chosen 
to avoid the Landau pole, where the Higgs mass is M# = 125 GeV and = 10 TeV with 
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FIG. 1: Evolution of 5 for Majorana neutrinos in the MSSM (upper plots) and in the UEDM 
(lower plots). The initial values 5 = ir/2, 5 = ir, and 5 = 3ir/2 are assumed, while the Majorana 
CP- violating phases p and a are marginalized. The values of # 12 , #13, #23 an d Am|i, in the 

la ranges from the global-fit analysis (for Am^ > 0) have been used as input 



121. 



R being the radius of the compactified extra dimension. Since the valid energy range in the 
UEDM is much smaller than that in the MSSM, the RG running does not develop as much. 
However, it should be noted that the RG running in UEDM is actually in the form of a 
power law, and thus can be more significant than in the SM and in the MSSM. 

It should also be noted that the Majorana CP-violating phases p and a have been 
marginalized over the range [0, it) in our numerical results. If the specific values of p and a 
are chosen, the variation of 5 will be even smaller. Therefore, we conclude that the leptonic 
Dirac CP-violating phase 5 is stable against radiative corrections in all the models under 
consideration, except for the MSSM with a large value of tan/3. In comparison, the Dirac 
CP-violating phase in the quark sector is stable even in the MSSM with a large value of 
tan /3, since the quark mass spectrum is strongly hierarchical. 



Now, we turn to the RG running 



:i.08±85?) n 



jehavior of 5 by taking the global-fit results 5 



ll| and 5 = (1.67lg$) tt 



121 ] as input. Since the present uncertainty is large, 
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FIG. 2: Allowed values of the leptonic CP-violating phase 5 (upper plots) and the Jarlskog invariant 
J (lower plots) for Majorana neutrinos at la C.L. with tan (3 = 10 (dark red or dark gray) and 
tan (3 = 30 (light red or gray) in the MSSM. The result of J in the MSSM with tan /3 = 50 is also 
given in the lower plots (yellow or lig ht g ray). The global-fit data from Ref. [TJJ are adopted for 
the left column, while that from Ref. 12| for the right column. 



we will choose the la range for illustration. In the upper plots of Fig. 2, the allowed regions 
of 5 at the superhigh-energy scale have been given in the MSSM. In the case of tan j3 = 30, 
one can observe that 5 is almost arbitrary within [0, 2ir) due to the large uncertainty of the 
input, so any predictions for 5 from a high-energy flavor model could be made consistent 
with the low-ener gy observations by the RG running. This is true for the global-fit results 
from both groups III, ll2j . In reality, any observable effects of CP violation should be related 
to the Jarlskog invariant J = Si2 c i2 s 23 c 23 s i3 c i3 s 5- Therefore, we also show the RG running 
of J in the MSSM for tan /3 = 10, 30, 50, in the lower plots of Fig. 2. It can be observed that 
J at a superhigh-energy scale could be quite different from that at the low-energy scale, in 
particular for tan/3 = 30 and tan/3 = 50. 
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III. RUNNING OF CP- VIOLATING PHASE: DIRAC NEUTRINOS 



The possibility for neutrinos to be Dirac particles has never been experimentally excluded. 
Moreover, it has been shown that the leptogenesis mechanism responsible for the matter- 
antimatter asymmetry in our Universe also works well in a different way for Dirac neutrinos 



2l| . Hence, we assume neutrinos to be Dirac particles, and give them masses through the 



sawa coupling 



22] 



coupling to the Higgs doublet — iiY v v^H + h.c. with Y v being the neutrino Yu 
matrix. It is convenient to write the RGE's of Dirac neutrino parameters as 

16tt 2 ^ = 2a u u + C vjL [(V/) oj + io{Y^)] , (9) 

where u = Y u Yj has been defined. The RGE's of k in the SM and the MSSM take the same 
form in Eq. (8), but with different coefficients a v and C ul , as given in Appendix B. Since the 
beta function for Dirac neutrino Yukawa couplings is currently not available in the UEDM, 
we consider only the SM and the MSSM. Similarly, as in the Majorana neutrino case, we 
find the RGE for the leptonic Dirac CP-violating phase 5 in the case of Dirac neutrinos 



S 23 C 23 S 13 S <5 



16vr 2 s 12 c 12 



£21 + ( C 12^32 S 12£3l) + 2 ^32 £31) 

S 13 



(10) 



where £y = (m 2 + m|) / (m? — m 2 ) has been defined. The relative error in the above equation 
is at the level of s 13 (Am| 1 /| Amfj-J) 2 ~ 10~ 4 . It is worth mentioning that the last term in 
Eq. (10) is comparable in magnitude to the second term, since the suppression by a factor 
of Amli/IAm^l is compensated by the enhancement from l/s 2 3 . Some general comments 
are in order: 

• The evolution of 5 is proportional to s s at all orders, so 5 will be kept unchanged by 
the RG running if s s = 0, namely, 5 = or 5 = ir. In other words, if leptonic CP 
violation is absent at low energies, it will never be generated by RG running. This is 
quite different from the Majorana case, where 5 can be radiatively generated via the 
non-vanishing Majorana CP-violating phases even if 5 = or 5 = it has been used as 
an initial condition. 

• Two qualitative differences between the SM and the MSSM should be noted. First, the 
tau- Yukawa coupling y\ = vr? T {\ + tan 2 /3)/w 2 in the MSSM is significantly enhanced 
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FIG. 3: Evolution of 5 for Dirac neutrinos in the MSSM for tan/3 = 10 (left plot) and tan/3 = 30 
(right plot). The initial values 5 = tt/2 and 5 = 3n/2 are assumed, and the values of 8 12 , #13, #23 
and Am| 1 , Am^ in the la ranges from the global-fit analysis (for Am^ > 0) have been used as 
input Q. 

for a large value of tan /3. Hence, the RG effect is more remarkable than that in the 
SM. Second, the coefficient C vl takes opposite signs in the SM and in the MSSM, 
indicating the evolution of 5 in opposite directions in these two models. 

To illustrate the RG running behavior of 5 in the Dirac neutrino case, we have shown 
in Fig. 3 two typical examples in the MSSM. In both examples, the initial values of 5 have 
been taken to be n/2 and 3tt/2, and the absolute neutrino mass is m 1 = 0.1 eV. The 
left plot is for tan/3 = 10, while the right for tan/3 = 30. Note that the beta function of 
5 is proportional to —s s in Eq. (9), where C ul = 1 in the MSSM. Therefore, 5 increases 
for 5 = 3ir/2, while it decreases for 5 = tt/2, as the energy scale evolves towards higher 
energies. This feature can be clearly observed in Fig. 3. Furthermore, the variation of 5 at 
any energy scale is quite small, compared to that in the case of Majorana neutrinos, where 
the arbitrary Majorana CP-violating phases play an important role in the evolution of S. 
As we have already mentioned, 5 will be kept unchanged if the initial values lead to s s = 0, 
so the trivial cases of 5 = and 5 = it have not been considered. 

Now, we continue with the global-fit results of S in Refs. 
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12J as initial values. The 



RG running of 5 in the MSSM for tan /3 = 10,30 and tan/3 = 50 have been shown in 
the upper and middle plots of Fig. 4, respectively. As before, the absolute neutrino mass 
m l = 0.1 eV is assumed. In the former case, the RG running effects are insignificant, which 
is in accordance with the results in Fig. 3. In the latter case, however, it is interesting to 
note that a wide range of values S G [0.27T, 1.87r] cannot be reached at the superhigh-energy 
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FIG. 4: Allowed values of the leptonic Dirac CP-violating phase 5 (upper and middle plots) and 
the Jarlskog invariant J (lower plots) for Dirac neutrinos at la C.L. with tan/3 = 10 (dark red 
or dark gray), tan/3 = 30 (light red or gray) and tan/3 = 50 (yellow or light gray) in the MSSM. 
The absolute neutrino mass m 1 =0.1 eV has been assumed. The global-fit data from Ref. [ill are 
adopted for the left column, while that from Ref. 12j for the right column. 



scale A = 10 14 GeV, no matter what initial value of 5 is chosen. The reason for this behavior 
is that the mixing angle 8 13 is approaching zero around A' = 10 s GeV. In the limit of an 
extremely small value of 8 l3 , Eq. (8) can be written as 



5 



16tt 2 



S 12 C 12 S 23 C 23 S <5 S 13 (^32 ^31 > 



(11) 
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where C vl — 1 has been chosen for the MSSM. Therefore, the RG running of S will be 
rapidly accelerated around A' = 10 8 GeV to the large- value region for s s < (i.e., 5 > n), 
while to the small- value region for s s > (i.e., 5 < it). This observation applies also to 
any initial value of 5. In fact, we have numerically checked the whole parameter region of 
5 G [0, 27r) at low energies, and found that only [0, 0.27r] and [1.87T, 27r) can be reached at 
high energies. However, the exact allowed range of 5 at high-energy scales really depends 
on the initial values of 5 and three mixing angles. For 5 = tt, the RG running of 5 will be 
absent, but 9 13 becomes negative above A' = 10 8 GeV, so we have to redefine 8 — > 5 ± tx 
to make 8 l3 positive, leading to 5 = or 27r at high-energy scales. In the lower plots of 
Fig. 4, the evolution of the Jarlskog invariant J is shown. Unlike the Dirac CP-violating 
phase 5 itself, the physical observable J evolves smoothly over the whole range of energy 
scales, as it should. For tan/3 = 50, the value of \ J\ can initially be as large as 2 %, it 
becomes vanishingly small at A = 10 14 GeV. One reason for this is that 5 shrinks into a 
small region around or 27r at the high-energy scale, as indicated in the middle plots of 
Fig. 4. Obviously, the evolution of the three mixing angles is also relevant here. 

IV. FURTHER DISCUSSIONS 

In Sees. II and III, we have examined the RG running behaviors of the leptonic Dirac 
CP-violating phase 5 in the cases of Majorana neutrinos and Dirac neutrinos, respectively. 
Now, we compare these two cases and summarize the main differences: 

• In the Majorana case, the two Majorana CP-violating phases are playing a crucial 
role in the RG running of 5. One can start from a CP-conserving scenario with 5 = 
or 7i at the low-energy scale, and end up with a CP-violating scenario even with 
5 = 7r/2 or 37r/2. In the Dirac case, the evolution of S is proportional to s s , so the CP 
conservation at the low-energy scale definitely implies that CP violation is absent at 
a superhigh-energy scale. 

• The mixing angle 6 13 could approach zero at some high-energy scale A' in both cases 
if a large value of tan/3 is assumed in the MSSM. On the other hand, there exist in 
the RGE's of 5 some terms inversely proportional to s 13 . Therefore, the RG running 
behavior of 5 will be dramatically changed around A'. Given the global- fit values of 
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S within the la range, it turns out that 5 could be arbitrary at the high-energy scale 
in the Majorana case due to the marginalization over p and a. In the Dirac case, 5 is 
found to be in two narrow ranges [0, 0.2tt] or [1.87T, 2ir) in the MSSM with tan /3 = 50. 

However, if a concrete mass model for Majorana neutrinos or Dirac neutrinos is assumed, 
the RG running of 5 may depend on the model details. In particular, when new particles 
or interactions come into play at some intermediate energy scale, the RGE's of the neutrino 



parameters are completely changed (23|. Hence, we have assumed that this is not the case 
in the previous discussions, at least below the cutoff scale. 

As we have mentioned before, many flavor symmetry models, which are intended for 
describing the observed leptonic mixing angles,predict the leptonic Dirac CP- violating phase 
5. For instance, it has been shown in Ref. 9] that 5 ~ 2n/3 (or An/3) and 5 ~ n/3 (or 
5ir/3) for different breaking patterns of the A 4 flavor symmetry in the type-I seesaw model, 
where three heavy right-handed neutrino singlets are introduced to realize the dimension-five 
Weinberg operator. If the vacuum alignment problem is further solved in the framework of 
supersymmetry, significant radiative corrections to these theoretical predictions of S could be 
possible. Thus, the leptonic Dirac CP-violating phase to be measured in neutrino oscillation 
experiments is related by the RG running to the theoretical prediction at the seesaw scale. 
On the other hand, the CP- violating and out-of-equilibrium decays of the heavy right-handed 
neutrinos can generate the lepton number asymmetry in the early Universe, which will be 
converted into the baryon number asymmetry via the SM sphaleron processes. In this 
case, the leptonic CP violation in neutrino oscillations can be associated with the matter- 
antimatter asymmetry in our Universe. 

V. SUMMARY 

Thanks to the recent measurements of 8 13 in the Daya Bay and RENO experiments, the 
discovery of CP violation in neutrino oscillation experiments seems to be promising if the 
leptonic CP violation really exists and the leptonic Dirac CP-violating phase 5 happens to 
be far away from or it. On the other hand, we have already had a preliminary result for the 
leptonic CP-violating phase 5 from the global-fit analysis of all kinds of neutrino oscillation 



experiments, namely 5 = (1.08to3i) tt [ll| and 5 = (1. 67^077-) ft [121 ] . Therefore, we are well 



motivated to study the RG running of 5 from the low-energy scale to a superhigh-energy 
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scale, where a unified model for fermion masses, flavor mixing, and CP violation is expected. 

In the case of Majorana neutrinos, we have introduced the dimension-five Weinberg op- 
erator to account for neutrino masses. The RGE of S has been derived analytically in great 
detail for the SM, the MSSM, and the UEDM, and a self-consistent approximation to it has 
been given as well. By a self-consistent approximation, we mean that the RGE of S has 
been expanded in terms of s\ 3 and Am^/IAm^l, and all the terms of the same order of 
magnitude should be preserved. It turns out that 5 is rather stable against radiative correc- 
tions in all these models, except for the case of a large tan (3 in the MSSM (e.g., tan f3 = 30 
together with a nearly degenerate neutrino mass spectrum). In this case, the Majorana 
CP-violating phases play an important role in the evolution of 5 such that a maximal phase 
5 = 7r/2 or 37r/2 can be radiatively generated at a superhigh-energy scale even if 5 = n (i.e., 
no CP-violating effects in neutrino oscillation experiments) at the low-energy scale. The 
evolution of § and the Jarlskog invariant J have been illustrated by taking the la global-fit 
results of 5 as input. 

In the case of Dirac neutrinos, we have derived the RGE of 5 in the SM and MSSM, 
and the self-consistent approximation to it has been made. Note that a nearly degenerate 
neutrino mass spectrum and the absolute neutrino mass m 1 = 0.1 eV are assumed in our 
analysis. The RG running effect of 5 can be neglected in the SM and in the MSSM with a 
small tan (3 (e.g., tan (3 < 10). However, 5 can be modified by more than 30 % for tan f3 = 30. 
The evolution of 5 and the Jarlskog invariant J have been examined by inputting the la 
global-fit results of 5. In the case of tan/3 = 50, 5 in the range of [0.27T, I.871"] is found 
to be unreachable at A = 10 14 GeV, since the mixing angle 6 l3 approaches zero at some 
intermediate scale (e.g., A' = 10 s GeV), which forces S to be in a large- value region for 
5 > 71 or a small- value region for 5 < it. At the same time, the Jarlskog invariant J becomes 
vanishingly small at a superhigh-energy scale. 

As we already know some information and will soon learn more about the leptonic Dirac 
CP-violating phase 5, it is thus meaningful to see how large it will be at a superhigh-energy 
scale. At such an energy scale, the leptonic Dirac CP-violating phase might be related to the 
quark Dirac CP-violating phase in a unified flavor model, or to the generation of matter- 
antimatter asymmetry in our Universe via the leptogenesis mechanism. In any case, the 
precise determination of 5 in the ongoing and upcoming neutrino oscillation experiments or 
at a future neutrino factory will shed light on the flavor dynamics at a high-energy scale. 
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Appendix A: RGE's for Majorana neutrinos 



1. The SM 



In the SM extended with the dimension-five Weinberg operator, the RGE for k has already 
been given in Eq. (3), while those for the Yukawa coupling matrices Y { of charged fermions 
(i.e., f = I for charged leptons, f = u for up- type quarks and f = d for down-type quarks) 
can be written as 



IQtt 



167T 



IQtt 
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dt 

dt 
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The relevant coefficients in Eqs. (3) and (Al) are Cf M = = Cfjf 
CSf = « = +3/2, and 
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with T 3M 



tr 



3 (y u y u f ) + 3 ( Y d Y l) + (Wi) ■ The RGE ' S for the su ^)c x SU(2) L x 



t/(l) Y gauge couplings g 3 , g 2 , and g l are given by 



167T' 



:^7i 
dt 



(A3) 



with (ftp 4 , 6| M , 6f M ) = (41/10, -19/6, -7). The quartic coupling A of the Higgs field appears 
in the RGE of k, which affects the evolution of absolute neutrino masses. It should satisfy 
the following RGE 
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It is worth mentioning that if the experimental uncertainties of the top quark mass M t and 
the strong coupling a s are taken into account, the SM vacuum could be stable up to the 



Planck scale A P1 = 1.2 x 10 GeV [24j, even for a Higgs mass M H = 125 GeV indicated by 



the recent results of the ATLAS and CMS experiments. 



2. The MSSM 



In the MSSM, the RGE's in Eqs. (3) and (Al) are still applicable, but the relevant flavor- 

i /-(MSSM _ r<^ 

1, , — C-„ 



universal coefficients are as follows: C^ 13 ™ = C^ 31 ^ 
C d M f 3M = 3, and 
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(A5) 



The RGE's for the gauge couplings are given in Eq. (A3), but with (bf 33M , bf ssu , bf SSM ) = 
(33/5, 1, —3) in the beta functions. As we can see from the RGE of k, the running neutrino 
parameters are determined by the charged-lepton Yukawa coupling matrix Y h especially the 
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tau-lepton Yukawa coupling y 2 . = m 2 {l + tan 2 (3)/v 2 , which could significantly be enhanced 
for a large value of tan /3. Such a unique feature can make the RG running of leptonic mixing 
parameters remarkable in the MSSM. 



The UEDM 



spacetime, so every 



25| . In the simplest 



In the UEDM, all the SM fields are promoted to a higher- dimensiona 
SM particle is accompanied by a tower of Kaluza-Klein (KK) modes 
UEDM with only one extra spatial dimension, which is compactified on an S 1 /Z 2 orbifold 
with radius R, the KK parity defined as (— l) n for the n-th KK mode is conserved after 
compactification. The mass scale of the first excited KK mode, i.e., /x = R~ l , has been 
constrained to be larger than about 300 GeV. 

If we extend the UEDM by an effective operator (£H) ■ k ■ (H T £ c )/2 to accommodate 
Majorana neutrino masses, just as in Eq. (1), then the effective Majorana neutrino mass 
matrix after electroweak symmetry breaking is M v = nv 2 with k = k/(jrR). The RGE of 
k now receives contributions from the KK modes, which are excited at the energy scale of 
interest. More explicitly, the RGE's for k and the Yukawa coupling matrices of the charged 



fermions are also given by Eqs. (3) and (Al), but with the following coefficients 25] 



a UEDM = «SM + s( _^ i2 _^ + a + 4T sm 
UEDM SM , „ / 33 2 15 2 orr sM 



a?" = «? M + s \--9t - -gi + 2T M 

a UEDM = a SM + s (_m g * _ Kfi _ + 2T s, 

UEDM ^SM , „ / 17 2 15 9 ^8 2 , rpSM I 

"d = "d + s I - ~ y 93 + M J ' ^ J 

and C*uedm _ C^ u (l+s) with "x" being any relevant subscript. Note that s = \p,j fi Q \ counts 
the number of excited KK modes for a given energy scale \x. In addition, the coefficients in 
the beta functions of gauge couplings turn out to be 

6 UEDM = 6 SM + ^ ^UEDM = 6 SM + 1 ? fe UEDM = fe SM _ * (A?) 

Finally, the RGE for the quartic Higgs coupling A is quite relevant in the UEDM, as in the 
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SM case. It has been found to be 25] 
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Appendix B: RGE's for Dirac neutrinos 

If the SM is extended with three right-handed neutrino singlets, then neutrinos acquire 
Dirac masses in the same way as the charged leptons and quarks do. At one-loop level, the 
RGE's of the fermion Yukawa coupling matrices read jzJ ] 
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with T^ M = tr 



3 (yX) + s to) + (YX) + fe f 



The RGE's of fermion Yukawa 



coupling matrices are the same as in Eq. (Bl) for the MSSM, but with different coefficients, 
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namely C f MSSM = +3 (for f = g) and +1 (for f ^ g), and 
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The RGE's of three gauge couplings g x , g 2 , and g 3 are the same as those in the case of 
Majorana neutrinos [see Eq. (A3)]. 
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